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Abstract 

Extended X-ray Absorption Fine Structure (EXAFS) experiments have been carried out on 
PrRu4Pi2 and PrOs4Pi2 to study the metal-to-insulator (MI) phase transition in PrRu4Pi2- No 
Pr displacement was observed across the MI transition temperature from the EXAFS data. In- 
stead, our EXAFS data clearly show that a Ru displacement is associated with this MI transition. 
The very high Debye temperature for the Ru-P bond (0£>=69O K) suggests that a slight rota- 
tion/displacement of relatively rigid RuPq octahedra leads to this small Ru displacement, which 
accompanies the MI transition at 62 K in PrRu4Pi2- 



The filled skutterudites have generated significant interest due to a broad range of novel 
physical properties, which include unconventional superconductivity, heavy fermion, and 
non- Fermi liquid behavior, as well as their potential thermoelectric applications in industry. 
Transport measurements show that most of these materials are metals with a very poor 
thermal conductivity, with severalphosphide skutterudites (such as CeFe 4 Pi 2 and CeRu 4 Pi 2 ) 
being semiconductors instead. jll Q Despite theoretical attempts to calculate the band 
structure of the filled skutterudites, a good understanding of the semiconducting behavior 
in these phosphides is still lacking. 

Recently, PrRu 4 Pi 2 was found to have a metal-to-insulator (MI) phase transition near 

\ n n 

Tmi of ~62 K. [3j Magnetic susceptibility and x-ray diffraction measurements |3J by Sekine 
et al. suggested that this MI transition is neither magnetic nor crystallographic in origin. 
Although valence fluctuations of the Pr ion might lead to this MI transition, a detailed X-ray 
absorption near edge structure (XANES) study of the Pr L/j-edge for PrRu 4 Pi 2 suggested 
that the Pr valence stays at 3 + above and below Tmi- Q Recent electron diffraction experi- 
ments by Lee et al. showed weak superlattice spots at 12 and 40 K, which were not observed 
at 70 K. 5] The observation was then attributed to a structural phase transition from Im3 
to another cubic phase, Pm3, as temperature decreases from 70 to 40 K. In the Pm3 phase, 
the Ru atoms are displaced along (1,1,1) directions. It was also suggested theoretically by 
Curnoe et al. that a symmetry lowering (1,1,-2) Ru displacement leading to a Pmmm phase 
may be responsible for this MI transition, though a (1,1,1) distortion mode leading to a 
Pm3 phase could not be excluded. JfJ Q] Band structure calculations also showed that a 
very tiny P distortion could open a band gap at the Fermi surface, leading to the MI phase 
transition in PrRu 4 Pi 2 . |s| 

Despite many experimental and theoretical studies of PrRu 4 Pi 2 indicating that a struc- 
tural distortion is responsible for the MI transition, there is still a lack of clear experimental 
evidence, which hinders a deeper understanding of this phenomenon. The main problem is 
that any structural distortion associated with the MI transition has to be very small, so that 
typical diffraction techniques cannot definitively observe it. Although recent synchrotron x- 
ray diffraction experiment on a single crystal PrRu 4 P 42 suggested a Pm3 phase at 10 K, 
|t| no high-T data were reported for comparison purposes. The extended x-ray absorption 
fine structure (EXAFS) technique, however, is ideally suited to detect small changes in local 
structure. We, therefore, carried out EXAFS experiments on both PrRu 4 P 42 and PrOs 4 Pi 2 , 
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because PrOs4Pi2 is a close structural relative of PrRu4Pi2, but is a metal below room 
temperature. Our data clearly show that a structural distortion of an unexpected nature 
accompanies the MI transition in PrRu4Pi2- 

EXAFS data at the Pr Ljjj-, Os Ljjj- and Ru K-edges were collected as a function 
of temperature at beam line 20-BM at the Advanced Photon Source (APS). Single crystal 
samples were ground into very fine powders and brushed onto scotch tape for the measure- 
ments (particle size ~ 10 //m). Energy-space data were reduced using standard procedures 







111 ] and the resulting fc-space data were Fourier transformed to r-space yielding peaks 



that correspond to different atomic shells. By fitting each r-space peak to a theoretical func- 



tion calculated using the FEFF8 program 



12j plus a Gaussian pair distribution function, 



we extracted the peak widths a and peak distances r. a provides quantitative information 
about the local distortions, including thermal vibration and static distortion, of a particular 
atomic shell. 

Fits to the Pr Lm-edge r-space data were performed for the first two atomic shells (Pr-P 
and Pr-Ru(Os)). Panel (a) of Fig. ^ shows the fit result to the 25 K data for PrRu4Pi2 
as an example of the fit quality. Since the filler atoms stay in a relatively large cage in 
the filled skutterudites, it is possible that the filler atoms might be displaced to an off- 
center position. Indeed, for PrOs4Sbi2, the larger value of a 2 at low temperatures indicated 



additional distortions which could be caused by a small Pr off-center displacement. 13] We, 
therefore, suspected that the MI transition in PrRu4Pi2 might be due to a Pr off-center 
displacement at low temperatures. However, a 2 for the Pr-P bonds, extracted from the 
fits (panel (b) of Fig. QJ) for PrRu4Pi2 and PrOs4Pi2 show no obvious anomalies as the 
temperature drops below 60 K. The rattling of the filler atoms is considered to be a localized 
vibration and, therefore, should be consistent with an Einstein model. 3| As shown in 
panel (b) of Fig. Q both sets of data can be very well fit to an Einstein model with no 
obvious deviation across Tmi- In fact, the static distortions, reduced masses and Einstein 
temperatures obtained in PrRu4Pi2 and PrOs4Pi2 are very similar, which suggests that the 
Pr local structure is essentially the same in both samples (static distortions ~0. 00003(5) 
A 2 ; reduced masses ~112(5) g/mol; Einstein temperatures ~130(3) K). Thus, we have no 
evidence for Pr off-center displacement in PrRu4Pi2 at any temperature, and consequently 
conclude that Pr displacement cannot be the origin of the MI transition in this material. 
We now compare the Ru K-edge data and the Os L in -edge data for the two materials. 
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Fig. |21plots the r-space data and the fit result at the lowest temperatures for each edge. The 
four nearest atomic shells (Ru(Os)-P, Ru(Os)-Pr, second Ru(Os)-P, and Ru(Os)-Ru(Os)) 
were used to fit the data. To account for a possible interference between further multi- 
scattering peaks and the Ru(Os)-Ru(Os) peak, two additional multi-scattering peaks at 
longer distances with small amplitudes were also included in the fit. 

Fig. El shows the temperature dependence of a 2 together with fits of cx 2 (T) to a Debye 
model for the Ru(Os)-P and Ru(Os)-Ru(Os) atom pairs. For the Ru-P and Os-P bonds, no 
anomaly in er 2 (T) was found across Tmi- The high Debye temperatures obtained for these 
two bonds (0 D = 690 K for Ru-P and 720 K for Os-P) suggest that the Ru-P and Os-P 
bonds are rather strong, making the RuPg and OsPg octahedra in PrRu4Pi2 and PrOs4Pi2 
rigid units. The most interesting feature in this figure is that the a 2 of the Ru-Ru pair starts 
to increase as T drops below 60 K (~Tm/)- If thermal vibrations dominated the broadening, 
a 2 should gradually decrease as T decreases, saturating at the zero-point-motion value. The 
a 2 for the Os-Os atom pair indeed follows this trend (O^ ~370 K). However, the increase in 
a 2 for the Ru-Ru pair at low temperatures clearly indicates that an extra distortion other 
than thermal vibration exists below 60 K in PrRu4Pi2- Fitting the Ru-Ru a 2 data from 
70 to 300 K yields ©£> of ~ 398 K with zero static distortion (the solid curve in Fig. OJ). 
Despite the excellent agreement with the T>70 K data, this fit yields a curve distinctly 
below the data for T<60 K. Attempts to fit the Ru-Ru data at all temperatures to a single 
Debye model plus a small static distortion (the dotted curve in Fig. EJ) yields a 0d of 
413 K and static distortion of 0.00018 A 2 . This model, however, fails to give a good fit to 
either the high- or low-temperature data. Thus there must be an extra Ru-Ru distortion at 
temperatures below 60 K. 

We now discuss the origin of this Ru-Ru distortion. In the undistorted lattice, all Ru-Ru 
distances are the same, so that any small expansion or contraction of the Ru framework 
that does not break this degeneracy cannot lead to a larger a 2 for the Ru-Ru peak. Thus, 
this extra distortion found for the Ru-Ru pair at low temperatures must come from a Ru 
displacement that breaks the degeneracy of the Ru-Ru distance. Lee et al. proposed that 
a Im3 to Pm3 phase transition is the origin of the MI transition in PrRu4Pi2- [5| Due to 
a (1,1,1) Ru distortion, the Ru sublattices then have two different sizes in the Pm3 phase. 
This splits the six nearest Ru-Ru atom pairs in the Im3 phase into 3 long and 3 short 
pairs, which could produce an extra broadening of the Ru-Ru peak in the EXAFS r-space 
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spectra. Although the low-T Pmmm phase suggested by Curnoe et al. changes the cubic 
Ru sublattice to tetragonal and also creates Ru-Ru distortion, ^,[?| recent inelastic neutron 
scattering measurements show a I\ ground state for Pr 3+ in PrRuiPi2, indicating that 
the assumption of a T 3 ground state in Curnoe's theory is not correct. Furthermore, the 
number of transitions of excited crystal field levels observed in the inelastic neutron scat- 
tering experiments below Tmi are consistent with the two inequivalent Pr sites in the Pm3 
phase. Therefore, the possibility of the low-T Pmmm phase can be excluded. 

To test the Pm3 model, we then fit the Ru-Ru peak in the EXAFS data with two Ru-Ru 
peaks. Note that EXAFS cannot normally resolve two peaks split by a small distance, as 
the fitting quality will be comparable for both a single-peak and a two-peak fit. However 
we can test whether the excess broadening observed at low T is consistent with the scenario 
suggested in the Pm3 phase under the assumption that the random static distortions are 
small. In this new fit, the number of neighbors was again fixed and the amplitude ratio of the 
two Ru-Ru peaks was set to be 1:1, each with degeneracy of three. Also, the a parameters 
for these two peaks were set to be the same, because the difference in the pair distances must 
be very small. Thus, only the Ru-Ru atom distances were allowed to change independently. 
As shown in panel (a) of Fig. 0J the low-temperature a 2 from the double-peak model are 
now lowered and very consistent with the 0£>=398 K Debye model which was used to fit the 
data from 70 to 300 K in the single peak model. The Ru-Ru peak distances obtained from 
the two fits are also plotted in Fig. 0] (panel (b)). The double-peak fit displays two Ru-Ru 
distances below Tmi, which are separated by about 0.022~0.030 A. This value is consistent 
with recently obtained low-temperature synchrotron x-ray diffraction results. [^] The result 
of the double-peak fit is clearly consistent with a Pm3 space group for T< T MI . Above T MI , 
the two Ru-Ru peaks collapse back together, confirming a high-temperature Im3 phase. 

As pointed out earlier, because the RuP@ octahedra are quite rigid units forming a corner- 
shared framework, a slight rotation, plus a tiny displacement, of the RuPg octahedra is 
needed to produce the observed Ru displacement. This inevitably leads to tiny correlated P 
displacements. However, because these P displacements are not principally along the radial 
directions of the Pr-P bonds, the subsequent Pr-P bond-length variations are too small to 
be observed in EXAFS (Panel (b) of Fig. EJ). 

Our EXAFS data clearly show that the MI transition in PrRu4Pi2 is associated with small 
correlated Ru and P displacements. Such atomic displacements may open a band gap at 
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the Fermi surface, leading to semiconducting (insulating) behavior at low temperatures. A 
band structure calculation by Harima et al. indeed showed that very small P displacements 
can introduce a band gap in PrRu 4 Pi 2 . 8] However, possible Ru displacements were not 
considered in that work. Thus, future theoretical calculations need to consider correlated 
displacements of Ru and P caused by a rotation/displacement of the RuPg octahedra. 

We also note that recent inelastic neutron scattering experiments of PrRu4Pi2 have been 
interpreted as evidence for a charge-density-wave produced by strong 4f-hybridization. 
Our results may help to illuminate this issue. Such a local-distortion-induced MI transition 
in PrRu4Pi2 may also be related to the semiconducting behavior observed in other phos- 
phide skutterudites, such as CeRu 4 Pi 2 and CeFe 4 Pi 2 . In this regard, the estimated static 
contribution to a 2 for the Ru-Ru atom pair from the Ru K-edge data for CeRu4Pi2 is about 
0.00093 A 2 at all temperatures; assuming a distortion of the Ru positions, this corre- 
sponds to a Ru displacement ~0.03 A (\J cr 2 tatic ) . A similar static distortion for the Fe-Fe 
atom pair in CeFe4Pi2 was also inferred. 

In summary, we present the first comprehensive local structure measurements of 
PrRu4Pi2- We find that a static Ru displacement is clearly associated with the MI transi- 
tion in PrRu4Pi2, while the P displacements must be associated with the motion of the rigid 
RuP6 octahedra. This Ru displacement found in PrRu4Pi2 suggests that the low tempera- 
ture phase in this material is likely to be Pm3. This EXAFS experiment clearly illustrates 
the importance of correlating the effects of local distortions with the transport properties of 
the filled skutterudites. 

This work was carried out under the auspices of the U.S. Department of Energy, Office 
of Science. The experiments were preformed at APS, which is funded by the DOE, Office 
of Sciences, Office of Basic Energy Sciences. Research at UCSD was supported by the US 
Department of Energy under Grant No. FG02-04ER46105. 
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FIG. 1: Panel (a): Pr L///-edge r-space data for PrRu4Pi2 at 25 K and the fit result. The Fourier 
transform range is from 3.0 to 10.5 A -1 , with 0.3 A -1 Gaussian broadening. The fit range is 
from 2.3 to 3.7 A. The high frequency curve inside the envelope is the real part of the Fourier 
Transform (FTr). The envelope is defined as: ±yj FT^ + FT 2 , where FT/ is the imaginary part 
of the Transform. There is a well-defined EXAFS phase shift for each peak, consequently the 
nearest Pr-P peak occurs at ~2.7 A and the Pr-Ru peak is at ~3.4 A. The inset displays Pr 
Ljjj-edge /c-space data at 25 K as an example to show the quality of the EXAFS data. Panel (b): 
the temperature dependent a 2 obtained from fits for the Pr-P bond in PrRu4Pi2 and PrOs4Pi2 
samples. The results of fitting a 2 vs. T to an Einstein model for the each sample are also plotted. 

FIG. 2: A plot of the Ru K- and Os L/jj-edges r-space data at 28 and 25 K along with the fit 
result for PrRu4Pi2 and PrOs4Pi2, respectively. The Fourier transform range is from 3.0 to 15.0 
A -1 , with 0.3 A -1 Gaussian broadening for both edges. The fit range is from 1.2 to 4.0 A for the 
Ru if -edge data; while it is 1.2 to 4.1 A for the Os Lj/j-edge data. The Ru-P peak is at ~1.8 A 
and the Ru-Ru peak is at ~3.8 A; while the Os-P peak is at ~2.0 A and the Os-Os peak is at ^3.9 
A. The insets in both panels plot the fc-space data for each edge, respectively, to show the quality 
of the EXAFS data. 
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FIG. 3: a 2 vs. temperature for the Ru(Os)-P, Ru(0s)-Ru(0s) atom pairs obtained from the Ru K- 
Os Luj-edges EXAFS data for PrRu4Pi2 and PrOs4Pi2, respectively. The Debye model fit to each 
set of data is also shown in the figure. For the Ru-Ru pair (solid circles), a 2 increases below 60K; 
here two Debye models are used: the dotted line is the fit to the data points at all temperatures 
and yields Bd of 413 K; the solid line is the fit result to the data points for temperatures between 
70 and 300 K, and gives @d of 398 K. The error bar for each data point, which is not shown, is 
smaller than the symbol size, except for one data point indicated in the figure. The error for the 
Debye temperature is ~20 K. 

FIG. 4: Panel (a): a 2 of the Ru-Ru atom pair obtained from the original single Ru-Ru peak fit 
and the double Ru-Ru peak fit as well as the 0£>=398 K Debye model. The inset shows the 
difference of these two sets of a 2 up to 100 K. Panel (b): The Ru-Ru peak distances obtained from 
the single-peak and double-peak fits. The error bars for all a 2 data and the peak distance in the 
single-peak fit are smaller than the symbol size and, therefore, not shown. 
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